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OF FINITE WIDTH MOVING THROUGH 
A GAS INITIALLY AT REST 
By Nathaniel B. Cohen 

SUMMARY 



The equations of motion and energy for the laminary houndary-layer 
flow in an expansion wave of finite width moving into undisturbed fluid, 
such as in a shock tube, were considered. Solutions in the form of 
infinite power series for velocity and local enthalpy functions were 
indicated, and the first three terms of each series were numerically 
evaluated. Validity of the numerical results was restricted to the 
region near the leading edge of the expansion wave. Skin friction and 
heat transfer were compared with values given by a solution which con- 
sidered the expansion wave as equivalent to a line discontinuity across 
which existence of isentropic expansion relations was assumed. These 
solutions were shown to be very different, qualitatively as well as 
quantitatively . Singularities in the flow field were discussed in 
regard to both the finite-width expansion-wave and the line-expansion- 
wave solutions. 


INTRODUCTION 


Extensive 'use of shock tubes for aerodynamic research has placed 
emphasis upon the effects of using a real, rather than a perfect, gas. 
The flow phenomena for an ideal fluid are easily derived and are given, 
for example, in references 1 to 3 . Deviations from the ideal flow may 
occur because the working fluid is imperfect and because fluid viscosity 
and thermal conductivity introduce effects of the shock-tube walls upon 
the flow. For moderate shock -pres sure ratios, imperfect gas effects may 
be neglected, but the wall effects may remain important; the present 
analysis considers only the fluid viscosity and conductivity. 

Deviations from ideal shock-tube flow are most easily seen experi- 
mentally through measurements of the shock -wave velocity. Experimental 
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attenuations in shock strength were investigated in references 1, 2, k, 

5, -and 6. Reference 6 also presented measurements of static -pres sure 
variations with time at a fixed position after passage of the shock wave. 
Experimental timewise density variations in the flow through the use of 
a chrono-interfercmeter were shown in reference 7- 

Various theoretical studies have been carried out as well. Refer- 
ences 1 and 4 considered a reduction in mass flow at the entropy discon- 
tinuity computed from the boundary-layer displacement thickness and the 
free stream corresponding to the unattenuated shock. Attenuation in the 
shock strength was then found by setting the mass flow through the shock 
wave equal to the reduced mass flow at the entropy discontinuity and com- 
puting the new shock strength. The unsteady laminary boundary layer in 
the hot gas behind the shock was derived in reference 2, and the resulting 
skin-friction and heat-transfer effects averaged across the assumed one- 
dimensional flow. The attenuation resulting from the waves generated at 
the entropy discontinuity by these boundary-layer effects were then 
computed. 

In reference 6, waves generated throughout the entire shock tube by 
skin-friction and heat-transfer effects were considered to affect the 
real gas flow, average skin-friction and heat -transfer coefficients being 
based upon incompressible, turbulent, steady-flow flat-plate boundary- 
layer values. The expansion wave was treated as a line discontinuity 
across which isentropic relations were assumed to be valid (zero -thickness 
expansion wave). Experimental data have shown excellent agreement with 
this theory. 

The assumption of the equivalence of steady- and unsteady -flow lami- 
nar boundary layers was shown to be rather inaccurate by the unsteady 
laminar boundary-layer theories of references 2, 8, and 9> which treated 
the exact boundary-layer equations for the flow behind a shock wave moving 
into an undisturbed fluid. Reference 9 extended the theory to laminar 
flow behind a zero-thickness expansion wave and presented solutions to 
the integrated equations of motion and energy for laminar and turbulent 
boundary layers behind both shock and zero -thickness expansion waves. 

The boundary-layer theories of references 8 and 9 have been incorporated 
into a theory predicting shock -wave attenuation in a shock tube (ref. 10). 

The turbulent boundary layer in an expansion wave of finite width 
(expansion fan) was considered on an equivalent steady-flow basis in ref- 
erence 11, but only recently has any attempt been made to investigate the 
unsteady boundary layer for this case. Reference 12 presents a solution 
to the complete shock-tube boundary-layer problem, including the fan, for 
laminar compressible flow. The integrated equations of motion and energy 
are considered and graphical and numerical isoclinal procedures are uti- 
lized. A feature of this theory is the use of the method of characteris- 
tics to obtain the solution, a necessity because singularities originating 
at the time of diaphragm burst propagate through the shock-tube flow field. 
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The present paper presents an approximate solution to the differen- 
tial houndary-layer equations of motion and energy for the unsteady 
laminar flow of a compressible fluid under the influence of the ideal- 
fluid expansion-fan pressure gradient. Results are compared with the 
results of the expansion-discontinuity analysis of reference 9 * Much 
of the material presented herein was submitted to the University of 
Virginia in partial fulfillment of the requirements for a Master of 
Aeronautical Engineering degree. 


SYMBOLS 


a 


C e 

C P 

c v 

F 

r n 

f 

f n 

Gn 

g 

Sn 

h 

h* 


ho 

k 

Npr 


local velocity of sound 
function defined by 

function defined by 

coefficient of specific heat at constant pressure 
coefficient of specific teat at constant volume 
particular function of nth-order differential equation of motion 
dimensionless stream function 

nth-order term in power-series representation of f 
particular function of nth-order differential energy equation 
dimensionless enthalpy function 

nth-order term in power-series representation of g 
local enthalpy 

enthalpy difference, h - hy- 
stagnation enthalpy 
thermal conductivity 
Prandtl number 


T i/ \^+ s / 

T e \°‘5/Ti + s\ 

h) sJ 
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p static pressure 

R gas constant 

q rate of heat transfer 


Qw - 


% 


2 re 
Pe a e /— 


<Jw = 

S 

T 

t, t 
u 




Pw N Pr /t5 


u. 


s 


X, X 


y 

7 

n 

V 


constant in Sutherland viscosity-temperature relation 

temperature 

time 

velocity of fluid in x- or x-direction 
wave velocity 
shock velocity 

velocity of fluid in Y- or y-direction 
coordinate along wall 


"incompressible" coordinate normal to wall, 

coordinate normal to wall 
ratio of specific heats, c p / c v 

v 

dimensionless coordinate, 

absolute viscosity 
kinematic viscosity 




0 P € 


7T- 
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| dj mens ionless coordinate, 1 — 

coordinate £ evaluated at trailing edge of expansion fan 
p density 

r shear stress 




i)r stream function 

0 ,cp dummy variable of integration 

Subscripts : 

n index of order of term in power serieB 

w quantity evaluated at wall (y = 0 ) 

a quantity evaluated in shock- tube region a, cold gas (see 

fig. 1 (a)) 

3 quantity evaluated in shock-tube region 3, hot gas (see 

fig. 1 (a)) 

s quantity evaluated in shock-tube undisturbed reference 

state e (see fig. 1 (a)) 

1 quantity evaluated in free stream (y -» 00) 

00 quantity evaluated in shock tube undisturbed low-pressure 

region » (see fig. 1(a)) 


Quantities written with a bar underneath are dimensionless with 
respect to the appropriate quantity in the undisturbed region e ; for 
u P 1 

example, u = — , p, = — , and so forth. 
a e — - L Pe 

Primes denote total differentiation of the functions f n (T]) and 
g n (q) with respect to the argument q. 
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THEORY 


The Prandtl boundary-layer equations are reduced to a system of 
two simultaneous, nonlinear, partial differential equations involving 
unsteady stream and enthalpy functions through a transformation of 
coordinates. The special case of flow in an expansion fan advancing 
into a fluid at rest is considered. Solutions in the form of infinite 
power series are assumed, and each partial differential equation is 
reduced to an infinite number of ordinary differential equations. 


Assumptions 

The assumptions necessary to the analysis which follows are: 

1. The Prandtl boundary-layer equations are assumed to be valid for 
the expansion fan, which is considered as centered, that is, originating 
at a single point on the x,t plane. A shock-tube wave diagram of the 
ideal-fluid flow for an expansion fan of moderate strength is shown in 
figure l(a). Figure 1(b) illustrates qualitatively the corresponding 
boundary layer. The point of origin of the fan at x = t = 0 (repre- 
senting, for example, the diaphragm burst in a shock tube) is a singular 
point and must be excluded from the analysis. The trailing edge of the 
expansion fan, separating the fan and the region of constant fluid prop- 
erties (region a, in fig. 1(a)), represents a discontinuity in the 
derivatives of the stream quantities and thus also violates the boundary- 
layer assumptions. Distortion of the trailing edge caused by varying the 
velocity and temperature within the boundary layer is neglected. 

In the region near the leading edge of 'the fan the velocity Is low 
and increases linearly in the ratio x/t from a value of zero at the 
leading edge. This situation appears analogous to the region near the 
forward stagnation point of a blunt body In steady flow and the Prandtl 
boundary -layer assumptions, through this analogy, are assumed to hold in 
this region. Therefore, the boundary-layer equations are assumed to be 
valid in the expansion fan except at the origin and along the free-stream 
trailing edge. 

One additional restriction on the use of the boundary-layer equa- 
tions must be noted for the special case of shock- tube flow. For a 
strong expansion process, part of the expansion fan occurs in the region 
where x > 0. The "no-slip" boundary condition at the wall requires that 
the temperature discontinuity separating the hot gas of region 3 from 
the cold gas of region a remain in the boundary layer extending back to 
the position of the diaphragm. This case is illustrated in the ideal- 
fluid wave diagram (fig. 2(a)) and in the corresponding boundary-layer 
formation (fig. 2(b)). At some time after the diaphragm burst, laminar 
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diffusion will have acted to smooth the discontinuity to some extent, "but 
large streamwise temperature gradients will still exist in the region 
near the position of the discontinuity in the free stream. This condition 
precludes application of the Prandtl "boundary-layer equations. Elsewhere 
in the region between x = 0 and the discontinuity, turning of the dis- 
continuity surface toward the wall will probably have reduced the stream- 
wise temperature gradient sufficiently for the boundary-layer equations to 
apply. Since the trailing edge of the fan always lies upstream of the 
discontinuity, it is assumed herein that streamwise temperature gradients 
large enough to invalidate the boundary-layer equations never exist inside 
the expansion fan. 

2. The coefficients of specific heat, Cp and c v , and the Prandtl 
number are assumed to be constant. 

3 . A linear viscosity-temperature relation is assumed. The constant 
of proportionality is defined so that the correct wall viscosity iB used; 
the error introduced elsewhere in the boundary layer is neglected. 

A. The wall temperature is assumed to be constant. Because walls 
are generally constructed of materials for which values of thermal con- 
ductivity and heat capacity per unit volume are high when compared with 
the values for the gas, this assumption, which has been shown to be 
valid for flow behind a shock in reference 9; should be reasonably good 
for the expansion-fan analysis. 


The Basic Equations of the Fluid Motion 

The differential equations governing the unsteady, two-dimensional 
flow of a compressible, viscous fluid under the Prandtl boundary-layer 
assumptions are the equation of continuity 


dp , d(pu) , d(pv) _ n 

St + "S -+ ~a^-° 


( 1 ) 


the equations of motion 
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the energy equation for constant Npj. and Cp 


8h 

8t 


+ 


u|ii + A = + 

ox 8y P \8t 




_JL_ 

pNpj. 8y^8yy 



( 3 ) 


and the equation of state 


P = 



oo 


The coordinate system is fixed to the wall, the origin of the x- and 
t-axes being at the origin of the expansion wave. ...The use of only two 
space dimensions, x and y, implies that interactions between adjacent 
wall boundary layers are neglected. 

Equation (2b) immediately yields the result that the pressure is 
constant across the boundary layer or that 


p(x,y,t) = p x (x,t) 


( 5 ) 


Equations (l), (2a), (3), and ( 4 ) then are to be solved for the four 
dependent variables u, v, p, h. The appropriate boundary conditions 
are 


u(x,0,t) = 0 
u(x,«,t) = u x (x,t) 


v(x,0,t) 

h(x,0,t) 


0 > 


M x 't) 


h(x, 00 y t) = h x (x,t)J 


( 6 ) 


In order that the equations may be more generally applicable, the 
enthalpy "at the wall is retained as a variable. 
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The linear viscosity -temperature relation discussed by Chapman and 
Rubesin (ref. 13) is employed; namely. 


iL = cv(x,t)JL (7) 

where the proportionality function (^(x, t) is obtained by requiring 
that equation (7) satisfy the Sutherland relation at the wall; thus. 


Cv(x,t) 


-5 fa* s\ 

X) x + s l 


( 8 ) 


A conventional transformation of coordinates and an unsteady stream 
faction are next introduced, the latter function identically satisfying 
the equation of continuity. The definitions are 


x = x 


Y ■ r y ay 

Jo P e 

t = t 


= 

P dy dY 


v= -v(^ + s/ o y ^^) -- 


> 


££(M. + 81 + 3 r\ 

P \dx dY dx dt/ 


(9) 


Equations (5 ) , (l ) , and the transformation and stream function of 
equations (9) are then simultaneously substituted into equations ( 2 a) 
and (3); the resulting equations are, respectively. 


d 2 ^ + d£_ d 2 ^ _ d£. dfV 7-1 _h_ ^1 + p l p l C w d% 

dYdt dxdY “ dx " 7 Pi dx Pg 2 Sy 5 


( 10 ) 
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3h,&^.Sh_S;l ! r > 3h_ 7-1 h /^1 , ^l\ ^1^1 ‘-'w c)^h , ^l^l^w/d^ilA 

3t oY dx dx SY 7 F[\dE SY ax / p^2 N Pr £y 2 Pg 2 \£y 2 / 

( 11 ) 


The appropriate boundary conditions on the stream function and enthalpy, 
written in x,Y,t coordinates, are 



N 


H) - u i(*^> 

dX /Y-*» 


i|r(x,0,t) = 0 
h(5c,0,t) = 3v(x,t) 

h(x,oo,t) = h-L(x,t) 


> 


( 12 ) 


Equations (10) and (ll) and the appropriate boundary conditions of 
equations (12) thus describe the boundary-layer flow in a -very general 
manner and require only the assumptions of constant Cp and Npp and 

a linear viscosity-temperature relation. Solution of a particular 
problem requires that the boundary conditions be given explicitly and 
that the differential equations be further reduced. The exact form of 
the boundary conditions will be shown to influence strongly the methods 
of simplification of the differential equations. 


Timewise Conical Similarity 

Similarity transformations ■ - Listed in appendix A are the potential- 
flow relations for flow in a centered expansion fan. The equations for 
Uj_ and hq are the free -stream boundary conditions for the boundary- 

layer differential equations and are functions only of the ratio x/t; 
thus, in a distance-time sense, the expansion-fan potential-flow variables 
are conical. If, in addition, the wall temperature is assumed to be con- 
stant, it is seen from equations (12) that all the wall boundary conditions 
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satisfy the coni cal requirement (that is, constancy), and it then appears 
reasonable to search for solutions containing this conical similarity. 

Similar stream an d enthalpy functions, f and g, respectively, are 
defined as 


*(x,Y,t) = U!(0^v e |€ f(|,n) 


(13) 


h(x,Y,t) = + h]_*(| )g(| ,tj) 


where 



(ifc) 


(15) 


The enthalpy at the wall hy is herein assumed to be constant. The 
velocity uy represents a wave velocity and is used so that, by employing 

the proper boundary conditions, equations ( 15 ) niay represent the flow 
behind a shock as well as in an expansion fan. For the flow behind the 
shock, the wave velocity is then u s j for the flow in the expansion fan, 

the wave velocity is considered to be the velocity of the leading edge of 
the expansion fan and is equal to -a e . The parameter | is the conical 
similarity parameter and is defined so that it is zero on the wave. The 
diffusion parameter r\ is closely related to that employed in references 2 
and 8. This relation may be seen by rewriting the parameter as 


1 (15a) 

‘ s) 


The length u^t - x is equal, to the distance from the wave, and the close 
resemblance of this parameter to that of the references is now evident. 
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Transformation of coordinates from x,Y,t to Z, ,r\,i through the 
use of relations (15) and simultaneous introduction of the stream and 
enthalpy functions (eqs. (13) and (l4)) reduce the equations of motion 
and energy (eqs. (10) and (ll), respectively) to the following: 


a u,c 3 l + [L . + + u 3f' 

^ v 6 3„5 3,2 [ 2t aAdT 2T 21 sr. 

f(5 - 1) + =L |£| fe, |£ 

[_ Hw \<^ 


. ^ a 2 ^ ) . 1 h * ^* 8 1 a & 




33w £1 as 


= 0 


(16) 


filHl' Pw 




^Pr S 


1 \6t) 2 > 




2£ 22^ d| 


— f + 
2& 


.)] + 

(S - 1) + =L §£ 

[d 

zl 

% 

d 


iil 




7-1 






= 0 


(17) 


Boundary conditions are transformed from eq.ua tlons (12) and are 


¥■) = 0 

Wt}=0 

fe) =1 

\ ' /T}-»oo 

f(S,0) = 0 

g(S, 0 ) = 0 

g(S,«) = 1 


( 18 ) 


The absence of the independent variable € from the differential equa- 
tions (lb) and (17) and from the boundary conditions of equations (18), 
except implicitly, indicates the existence of solutions exhibiting simi- 
larity in | and tj as was assumed in equations (13) and (l4). The 
number of independent variables has therefore been reduced to two. 
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Application to the centered expansion fan .- The special case con- 
cerning flow in the boundary layer behind a shock wave advancing into a 
fluid at rest leads to a result identical to those of references 2 and 
8. This case is discussed in appendix B. The differential equations 
governing the boundary -layer flow in the centered expansion fan are 
found by substituting the free -stream relations of appendix A (.equa- 
tions (Al) to (A5) and (AT)) into equations (l6) and (17) with the 
further qualification that 


u w a e 


I 



(19) 


Results of this substitution are 

3y.pi 


a5f(a + 5_ijf + 

3„ 2 V 2 27 + 1 


i(i .mlA 7 ' 1 2l + 

«€\ r + i/ a n 3 V r+1 / 

jA* I) + [ (l - 11 - Tfi 5 ^ + s fe)f + 






g = 0 


cA r+i 


5Z=i /_ 




8 2 g 

M Pr Sti 2 


< 7 - + ( L ■ Ttl 5 ) [ (1 ■ ^ • 2 ^ + 

|sfa + 1 - 2 - r 4f + 


Z_zA\ |2 
7 + 1/ 

7 “ 


^ 2 &] + fit? 


( 20 ) 


(ft - U - 


2 t 8f 


8^\2 T 2 7TT" T TTT 8F; ‘ " 7 + r^J 

~ I | (1 _ Z_I_i|\ 2 (| - 1) £ — ||£lg + 2 y ~ A (i - 1) - 

+ 1 S \ 7 + l 5 /[_ 7 + 18 t]J 7 + 1 


V - 


2 .df 


7 + 1 8 t] 


{i - ^ + AtrA}} - 0 


( 21 ) 
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where 


S*. _ &T' 5 /v l 1\ 

C e \T„ + S j 


The boundary conditions upon the functions f and g are the same as 
those given by equations ( 18 ). 

Because of the complexity of these partial differential equations, 
solutions are sought in the form of infinite power series, which are 
written in powers of the conical coordinate | with the coefficients 
as functions of q only. Thus, the solutions are assumed as 


= XI £%(n) 

n=0 


g(l^T)) = | n gn(Tl) 


( 22 ) 


Because both of these partial differential equations possess singular 
points in the region £ j? 1 (see section entitled "Results and 
Discussion"), the infinite series can only represent the solutions in 
the region g < 1. The first few terms of each series are considered to 
approximate satisfactorily the exact solution in the region near the 
leading edge £ « 1. 

One further assumption is desirable. The enthalpy at the wall has 
already been assumed to be constant but arbitrary. The most common 
shock-tube physical arrangement dictates that the wall temperature before 
passage of the wave must be equal to the temperature of the undisturbed 
fluid. Thus, if the wall temperature is assumed to remain unchanged after 
passage of the wave because of the high thermal conductivity and heat 
capacity of metal walls, then 


^w = h e 



( 23 ) 



NACA TN 3943 


15 


Substitution of equations (22) and ( 23 ) into the equations of motion 
and energy (eqs. (20) and (21)), application of the boundary conditions 
(eqs. (l8)), and equation of coefficients of like powers of | succes- 
sively leads to an infinite set of ordinary linear differential equations 
and boun dar y conditions for each function, the first three of which are 



f L w 


f 


tn 


2 


+ ^o" 

- f 0 * = “I 

(24a) 

+ 

- 2f 1 * = Fi(ti) 

(24b) 

+ 2 f 2 

- 3f 2 * = 

(24c) 


f 0 (0) = f Q '(0) = 0 

f 0 '(“) = 1 
f n (°) = V<°) = f n'(“) = 0 


(25a) 

(n ^ l) ( 25 b) 


fr - + ^0 ' - «o - - 1 


N- 


Pr 


§; + ‘ 

5^ + a'fe' • 5e 2 = 
go(o) = 0 
So(“) = 1 

* 


( 26 a) 


( 26 b) 


( 26 c) 


( 27 a) 


(ntl) 


Sn(°) = Sn(°°) = 0 


( 27 b) 
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with the following definitions : 


F l^) j + 7 l f o" “2 7 + l f 0 f 0 f o' + y + i( f 0') 2 +2 7 + 1°0 "7 + 1 


(28a) 


F 2 (ti) = 


- ^ 1 ' - 5 7 fr f0,f1 ' - 1 7TT f ^" - 1 - 

+ + {5riJ ( 8 ° ‘ x ) (28t) 


OxtuK- ,r-5 y + -^-(f 0 ") 2 - 2-i-Joso' + 

1MA 2 7 + 1 2(7 + 1) 2®° 7 + l k ^ 27+I U&U 


7 

2 


— f o’ g o “ 7T7 f ° 


(28c) 


G 2 (n) = 


_i 7 _ , 7 ....t-.. 7 -..g 1 . 1 - 2 — fo gl > - f lg0 ' + 

7 + 1 w Pr 2(7 + 1) 1 27+1 01 27 + 1 1 

2 _^_f -g + — i— -f-L'go §— f^ + 2 _£_.f 0 ,, f 1 " - 

7 + 1 u x 7 + 1 x u 7 + 1 x 7+1 x 


3 T 2 + 187 “ 13 So" , T 2 - 147 + 13 T]„ , c w/ - , 

+ — ^ 0 - oUTTT)trTT) I 0 


3/ 2 \/7 - 1 


4(7 + 1)2 4(7 + 1)2 2 2\7 + ly\7 + \ 

1 2 37 + 1/4. .1,2 , 3 2 7 - I4. 

5 — Tn (f °) + 2 7nrn foB ° 


(28d) 


The primes on a symbol indicate differentiation with respect to the 
argument T) . 

Numerical integration of the preceding equations can he accomplished 
with no special difficulty. The momentum, and energy equations for the 
nth functions are decoupled and hence may be solved independently. The 
particular integrals F n (ri) and G n ( T| ) are dependent only on previous 

solutions. 
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Equations (24) to (28) are solved for f n (ri) and g n (Ti) and their 

respective derivatives , for n = 0, 1, and 2. These solutions then 
make possible approximation of the stream and enthalpy functions by the 
following relations which are valid in the vicinity of the leading edge: 


f(i,ri) « f 0 (^) + 6fi(n) + s 2 f 2 (n) (29a) 

g(s,ii) « go( T i) + Ssi(n) + 5 2 g 2 (n) (29b) 

Velocity and enthalpy profiles are given by 

f 0 ’(Tl) + Ifi’(n) + S^’Ol) (30a) 

= g(s,n) ~ eo(n) + lgq(n) + l 2 g 2 (fi) (30b) 

V h i - K 


Skin Friction and Heat Transfer 
The skin friction is given by 

Tv = (^) y=0 ‘ 


(31) 


This relation is transformed to the g,T),t coordinate system through 
the transformation equations (9), (13)> (l4) , and (15); the transforma- 
tion yields the dimensionless friction 


T w = 


•w 



£l£l c - 


w 


% 

ft 




”* TJ=0 


(32) 


Finally, for the special case of expansion-fan boundary layer, from equa- 
tions (23), it is seen that 
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where 



Sk 

ft 






Jt ]=0 


(33) 


a 2 f(l ' T l ) = f 0 "(0) + |f 1 "(0) + l s f 2 "(0) (35a) 

. 3,12 Ji=o 


nip heat transferred per unit time per unit surface area from the 
wall to the fluid q^ is given hy the relation 


^ - -"“(ll (54) 

The transformation to the (■ jt^t system is accomplished in the same 
manner as for the skin friction and yields the dimensionless heat 
transfer 


^w 



■£l£l. C w _ 1 _ 


N. 


Pr 


^L* 

Sg(l^) 

7-1 



-*T]=0 


(35) 


The expansion-fan wall heat transfer is obtained from equation ( 35 ) and 
is 


where 


^w 



«!P 

ro 

H 

dg(6,Tl) 

/f 7 - 1 

8t] 


a 


dg(S jTi) 


8 t] 




=0 


go'Co) + Up'Co) + l 2 g 2 '(o) 


(36) 


( 36 a) 
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Transformation to the Physical Plane 

The transformation to the physical coordinate system is carried out 
by using the profile functions. The incompressible normal coordinate Y 
was defined by equations (9). The reverse transformation is 


y = 


0 p(x,Y,t) 


-ST— dX = 


r 5 

P x (x,t) Jo 


h(x,Y,t) 

hc- 


dY 


( 37 ) 


Since the local enthalpy is given as a function of | and T] by 
equation (30b), the integration of equation (37) is most easily accom- 
plished in those coordinates. The transformation equations defining | 
and t} (eqs. (15)) are used to define the differential dY as 


dY = (tor + (58) 

The integration (eq. ( 37 )) is to be carried out at constant it and t 
and thus at constant | and tj as a result of combining equations (38) 
and (37), the physical normal coordinate y is 


Y = 



an 


( 39 ) 


Substitution of the enthalpy function from equation ( 30b ) , rearrangement 
of the terms, replacement of t by its equivalent t, and integration of 
equation (39) yields for = h £ : 


h * 

1 —1 

P, P 


)jv e £t E-i 5.1 


n T 1 p TJ p TJ 

g^^dT] + s J gqCr^dq + £ 2 J 0 Sg^)^ 


(bo) 


The quantity ■ is the physical compressible counterpart of the 

incompressible similarity parameter tj and is a function only of | 
and T). 



20 


MCA IN 3943 


RESULTS AND DISCUSSION 


Numerical Results 

Each of the homogeneous differential equations corresponding to the 
complete equations (24) and (2 6) was transformed to Hermite's differential 
equation and complementary solutions in the form of Hermite polynomials 
were obtained. (See ref. 14.) Except for the n = 0 case, particular 
solutions appeared to be unobtainable in closed form. Rather than combine 
the Hermite polynomials with particular integrals found by numerical pro- 
cedures, it was decided to treat by numerical methods the complete equa- 
tions for the cases n ^ 1. 

In reference 15, the laminar boundary layer on a flat plate uniformly 
accelerated from rest was investigated. The resulting momentum equation 
and boundary conditions for the initial motion of the plate are equivalent 
to those of the present analysis for n = 0 (eqs. (24a) and (25a)) and the 
complete solution (from ref. 15) is therefore 



It can be shown that the first term on the right-hand side of equation (42) 
is the particular integral and the two remaining terms are the exponential 
form of the appropriate Hermite polynomial, the complementary solution. 


The energy differential equation is of the same form as the momentum 
equation, and the following solution may be deduced: 


Bq(ti) = 


-N 


Pr 




Th/%^ -^(11/2) 2 Npi>n 2 + 2 


N p tj/2 


f 


ft 


de 


( 43 ) 


These functions satisfy the respective boundary conditions and are given 
in table i for N Pr equal to O.72. 
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Numerical integration of the higher order differential equations was 
carried out on the Bell Telephone Laboratories X-66744- relay computer at 
the Langley Laboratory by using a modification of the Runge-Kutta fourth- 
order method given in reference l 6 . (The ratio of specific heats 7 was 
set equal to l.lt-.) Results of these computations are listed in tables II 
and III. The functions most important to the analysis are fn'C^) and 
g n (Ti), the nth order velocity and enthalpy profiles , respectively; plots 
of these functions are given in figure 3 - 

The actual velocity and enthalpy profiles in the |,Tj,t coordinate 
system are approximated by equations ( 30 a,) and ( 30 b) and have been com- 
puted for values of | equal to 0, 0.1, 0.2, 0 . 3 , and 0 . 5 . The results 
of this computation are shown in figure 4. Since the power series for 
'both velocity and enthalpy profiles neglect terms of order |3, the results- 
are considered to be quantitatively correct only in the range 0 ^ 0 . 3 .' 

Extension of the results up to | = 0.5 is included for qualitative indi- 
cation of trends. Skin-friction and heat-transfer functions are shown in 
figure 5 'with three approximate formulas for each, corresponding to taking 
one, two, or three terms of the power series. 


The transformation integrals were computed by graphical integration 
of the appropriate functions given in tables I to III, and are listed In 
table TV. The function defined by equation (40) is computed for values 
of {■ of 0, 0.1, 0.2, 0 . 3 , and 0 . 5 , and the results are plotted In fig- 
ure 6 . Theoretical velocity and enthalpy difference profiles in the 
physical x,y,t coordinate system have been computed and are given in 
figure 7 for values of | of 0.1, 0.2, 0 . 3 , and 0.5 with the Infinite 


series approximated by the three term serie s . T he normal coordinate 
parameter yj\jv e t Is used In place of yjjv^t so that adjacent pro- 
files may be more easily distinguished. 


Discussion of Results 

As mentioned previously, the zero -order momentum equation for the 
present case was found to be Identical to that of reference 15 , the ini- 
tial motion of a flat plate uniformly accelerated from rest. This result 
is not surprising in view of the fact that the present case may be con- 
sidered to be a form of constant acceleration; that is, the potential- 
flow velocity varies linearly with the conical coordinate | . The 
respective coordinates of the reference are also closely related to the 
coordinates | and tj used herein. 

Figures 4 and 7 are apparently contradictory; figures 4(a) and 4(b) 
show profiles becoming increasingly full as £ is increased whereas the 
corresponding plots of figures 7 (a) and 7 (b) indicate velocity and 
enthalpy difference profiles behaving in the opposite manner. This 
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condition is the effect of distortion of the normal coordinate by plotting 
it in the incompressible T} -plane as opposed. to plotting it in the physical 

y jf 7 t -plane, and no contradiction actually exists. 

For | equal to 0 . 5 , both velocity and enthalpy difference profiles 
in figures 4 and 7 show a slight tendency to oscillate once near the free- 
stream boundary, the latter profile actually exceeding a value of unity 
inside the boundary layer. The trend indicates that the velocity also 
will exceed unity within the boundary layer for values of | larger 
than O.5. Because the amplitude of oscillation is extremely small and 
the estimated error of the present solution for the three -term-series 
approximation is of the order of these ^results are open to question. 

However, velocity profiles exceeding unity have been encountered in 
steady-flow theories under conditions of a strong favorable stream pres- 
sure gradient and high ratio of wall temperature to stream temperature _ 
(ref. 17); such stream and wall conditions are present near the leading 
edge of the expansion fan. Thus, although these profiles must be ques- 
tioned in view of the limits of validity of the present series solutions, 
the velocity profile, at least, is typical of. certain steady-flow velocity 
profiles in fluid undergoing similar processes. 

The heat transfer given by the three-term-series approximation 
(fig. 5(b)) apparently reaches a maximum at ? » 0.4 and then decreases. 
This result may be qualitatively explained by considering the difference 
between the wall enthalpy and the stream stagnation enthalpy, which 
forms a sort of heat-transfer potential. The stream stagnation enthalpy 
h 0 is derived from the expansion-fan relations of appendix A and is 

given by - 


ito,l 



1 + 2 LrL 4 fe 2 - 25) 

7 + 1 \ / 


( 44 ) 


The dimensionless potential, h^ - q for the ease h^. = h £ , is 

plotted in figure 8. The rate of heat transfer follows the same general 
trend as the total enthalpy difference for an increase in the conical 
coordinate near the leading edge, 


If this line of reasoning is followed, one might expect that, in the 
region where the stagnation enthalpy difference becomes very strongly 
negative, the rate - ' of heat transfer might also become negative and would 
thus indicate heat transfer from the fluid to the wall. This is indeed 
the case with the integral solution of reference 12 where the heat- 
transfer function goes from positive to negative at a value of the conical 
coordinate | 1.5. The present solution cannot be carried to that 
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extent, but examination of the enthalpy function of the second power of 
£ (S2( T l) fig- 3(b)) indicates a reversal of sign in the derivatives 

at the wall (rj = 0); thus, a tendency for the heat transfer to reverse 
is shown. This term in the power series increases in importance with 
increasing £ and lends support to the belief that regions of negative 
rates of heat transfer (from fluid to wall) exist beyond the region of 
validity of the present solution. 


Comparison With Solution of Zero -Thickness Expansion Wave 

Reference 9 considers the expansion fan as a wave of zero thickness 
propagating into fluid at rest at a velocity equal to -a e . The wave 

diagram for the special case of shock-tuDe flow under this assumption is 
shown in figure 9(a); the true fan is illustrated in figure 9(b). The 
zero thickness wave coincides with the line g = 0. Isentropic expansion- 
fan relations are assumed to hold across the zero-thickness wave, and the 
entire region from the wave to the temperature discontinuity is assumed to 
be a region of constant stream properties and corresponds to region a in 
the true shock-tube case. Comparison of the results of the present analy- 
sis wit h those of reference 9 is carried out on the basis of equal poten- 
tial flows. That is, if the value of g at the trailing edge of the 
expansion fan is given, free-stream variables in the constant poten- 

tial flow region a, (that is p , u , and so forth) may be found from 
appendix A (eqs. (Al) to (A5) evaluated at g te ). These variables are 

also the free-stream variables for the corresponding zero -thickness 
expansion-wave solution. Since the present analysis does not consider the 
region g > g^ (region a), comparison must be restricted to the region 

0 = i = 6-fce for each case, although the results of reference 9 are 

extended Deyond this region in figure 10. The asterisks appearing on the 
curves represent the appropriate positions of g^ e . 

Skin-friction and heat-transfer functions are shown in figure 10 for 
it e equal to 0.1, 0.2, 0.3, and 0 . 5 , the present solution corresponding 

to the three-term (quadratic ) power-series approximations (eqs. (33a) and 
(36a) ) . It is apparent that the two solutions are very different; the 
skin friction and heat transfer for each zero-thickness expansion solution 
approach infinity at the wave 1=0, whereas the corresponding functions 
for the expansion fan vanish. One might consider this behavior analogous 
with that of a flat plate and of a blunt body in steady flow. Skin fric- 
tion and heat transfer at the leading edge of a flat plate are infinite 
whereas, at the forward stagnation point of a blunt body, friction vanishes 
and heat transfer is finite, depending upon wall temperature . Whether 
these discrepancies are Important in computing shock attenuation, however, 
will depend upon the strength of the wave in question. For weak waves, 
although the details of the flow within the narrow fan have been shown to 
be very different, effects upon shock attenuation may be small as they 
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represent integrated wall shear and heat transfer far downstream. Here 
the steady-flow analogy is with the flow well downstream of the leading 
edges of' the flat plate and the blunt body where the leading edge is not 
important. For strong waves , however , the flow within the fan may 
seriously influence the shock wave. It should be pointed out that 
although the computations of reference 9 were carried out for strong as 
well as weak expansions, the limitations of“the analysis were therein 
recognized. 

It is of interest to note that, at | = in the zero-thickness- 

wave case, the skin-friction and heat-transfer functions are approximately 
one -half the corresponding expansion-fan values. This relation is more 
evident if skin-friction and heat-transfer coefficients defined by 



(45) 


(46) 


are considered. For the present solution, combination of equations (35) 
and (45) and of ( 36 ) and (46) yields, respectively, 


T w 



(47) 


%r 



(48) 


The corresponding coefficient form from reference 9 is independent of 
coordinate % and is a function only of expansion strength |^ e . These 

coefficients are shown in figure 11; the values from reference 9 are 
approximately one-half those of the present solution. This relation is 
exact for the limit |^- e = 0 (an infinitesimally weak wave) where the 

following values apply: 
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Present solution 

Reference 9 

T w 

2fi 

£ 

% 

2 *Er 

1~ 

* 


The zero -thickne s s -wave case is discussed, further in appendix B. A 
more detailed comparison of the finite-width and zero-thickness expansion 
waves is made in reference 12 . 


Singularities in the Flow Field 

The equations of momentum and energy were transformed to the |,t|,€ 
coordinate system (eqs . (l 6 ) and ( 17 )) and then specialized for the case 
of the centered expansion-fan potential flow (eqs. ( 20 ) and ( 21 )). It is 
of interest to note that the aforementioned equations indicate mathemati- 
cal singularities in the flow field. This condition can be observed in the 
expression 


(t - 1 ) - ^ §£ 



(5-1) 

L 3 ii 


7 + 1 ch} 

_ 


( 49 ) 


which is the coefficient of the highest order derivative with respect to 
| in each of the partial differential equations. For 1 5 | =1+Hi> 

this term vanishes somewhere in the boundary layer and, as a result, the 
solution became nonanalytic (that is, a singular point). In the x,y,t 
coordinate system, the singularity occurs at a position 



'singularity 


= Hq || = 3 i(x,y,t) 


It is evident, since u is always positive (in the positive x-direction), 
that the singularity always appears in the region x ^ 0 and is actually 

within the limits 0 ^ - 2 — ^ u. . 

a e t -1 

The appearance of these singularities is discussed in reference 18 
in regard to the impulsive motion of a semi-infinite flat plate in an 
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incompressible fluid. David Adamson has made an independent investigation 
of this case at the Langley Laboratory and he argues that, although the 
vorticity generated at the plate surface at the start of motion diffuses 
normal to the surface at an infinite rate, the Prandtl boundary-layer 
assumptions allow no diffusion of vorticity in the streamwise direction. 

At the origin of motion, the flow field of this case consists of the 
undisturbed fluid ahead' of the leading edge, and the fluid behind the 
leading edge through which the vorticity of starting has already spread 
in a direction normal to the surface. A line singularity separates these 
regions and, as time progresses, this singularity is swept downstream. 

The singularities in the expansion fan may be considered in an analo- 
gous fashion as originating at the point x = t = 0 (the diaphragm sta- 
tion at burst in shock-tube flow, for example) with vorticity generated at 
this point and diffused immediately only in the normal y-direction. As 
time progresses, the singularity is swept downstream (in the positive 
x -direction) . When the expansion fan is supersonic (£te > T), these 

singularities appear in the fan itself; thus, the present series solu- 
tion is limited as already discussed. 

The singularities exist, however, even when the expansion fan is sub- 
sonic, that is, 5 te < 1. The condition of a nonanalytic solution must 

still occur in the region where 5=1 even though this region may be down 
stream of the expansion-fan trailing edge bercause the physical reasoning 
describing the diffusion of vorticity is qualitatively independent of 
expansion-fan strength. As a result, one would expect the analysis based 
upon a zero-thickness expansion wave to exhibit singular behavior in the 
region 5=1. The fact that no such behavior was evident in the data of 
reference 9 is attributed to the fact that the origin of the wave (for 
example, the diaphragm burst) was neglected therein when the flow was 
assumed to be steady with respect to a coordinate system fixed to the wave. 
The zero-thickness wave, as applied to the shock tube in references 6 and 
10, therefore implied a wave which originated at t = -» , propogated at 
a constant velocity equal to -a e , and passed through the point x = t = 0 

on the x,t diagram (fig. 12). Consideration of the zero-thickness expan- 
sion wave by the present method (see appendix B) does result in the 
expected mathematical singularities. 


CONCLUDING REMARKS 


A solution to the laminar boundary-layer equations that is valid near 
the leading edge of an expansion wave of finite width has been obtained. 
Application of the results to the calculation of real gas flow in a shock 
tube by the method of NACA Technical Note 3375 involves the use of the 
skin-friction and heat-transfer results along with like functions appro- 
priate to the cold-gas constant-potential-flow region. The present 
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solution could also be adapted for application to the method of NACA 
Technical Note 3278 hy computation of the velocity normal to the wall. 

For this application, solutions in the region between the fan trailing 
edge and the temperature discontinuity, including the associated singu- 
larities, also would be required. 

The zero-thickness expansion-wave solution of NACA Technical 
Note 3712 indicates the skin-friction and heat-transfer results to be 
very different from those obtained herein. In principle, then, substi- 
tution of a zero -thickness wave for a true expansion fan in computing the 
boundary layer must be questioned both on grounds of accuracy near the 
leading edge and of omission from the consideration of the singularities 
in the flow field in the region downstream of the position of the origin 
of the wave. In application to the shock-tube case, however, the results 
of NACA Technical Note 3278 indicate that the hot gas compressed by the 
shock is the primary factor influencing shock-wave attenuation for moder- 
ate pressure ratios; thus, in practice, assumption of a zero -thickness 
expansion wave may not seriously affect the theoretical attenuation 
results . 

No experimental boundary-layer data for laminar flow in a centered 
expansion fan are available. Local heat-transfer rates within the fan can 
be measured, and such measurements may substantiate the present theory. 
Data on shock -wave attenuation are plentiful, but little information 
regarding the fan may be deduced since the effects of the shock-compressed 
gas generally predominate and flow Reynolds numbers are frequently in the 
turbulent range. Future application of the integral solution to turbulent 
boundary layers should afford both direct and indirect checks on that 
theory and should lead to more definite conclusions regarding the present 
theory . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 3; 1956. 
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APPENDIX A 

FREE-STREAM RELATIONS IN THE EXPANSION FAN 


The equations defining the flow parameters of a fluid undergo ing 
acceleration in a centered expansion fan are listed, for example, in 
reference 1 to 3 and are reproduced below in terms of the coordinate 


| = - 2 _ + 1 : 

a_t 


% = — - 2 -rS 
1 7+1 

(Al) 

Sq = 1 - 7 

7+1 

(A2) 

CM 

II 

(A3) 

2 7 


El * {hf* 

(A4) 

2 


£1 ■ (%) r " 1 

(A5) 


The stream viscosity is written by using the Sutherland relation as 


„ - n - * a) 

-1 M’g \T e j \T-j_ + S) 


(A 6) 


Equation (a 6) is combined with the equation of state (eq. (4)) and the 
definition of C v (eq. (8)); from the resulting relation are obtained 
the viscosity terms in the differential equations of motion and energy 
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vhere 


x 

c e 


/ T l\°‘ 5 /' T e + s \ 

\Te ) \T]_ + SJ 


Equations (8) and (A8) may be combined to yield 


Oj, /l„\ 0 -5/ I 6 + s\ 

c e \^J V>V + £ V 

For the special case = T e , it is seen that 



(AS) 


(A9) 


(aio) 
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APPENDIX B 

BOUNDARY LAYER BEHIND SHOCK OR ZERO -THICKNESS 
EXPANSION WAVE MOVING INTO FLUID AT REST 


The flow behind a shock wave or a zero-thickness expansion wave may 
be considered by utilizing the theory of this paper modified for constant 
potential flow. The x,t diagram of figure 9(a) shows both cases. 

Although the shock wave advances into the low-pressure region «>, region e 
remains the reference for both solutions . The wave velocity u^ is equal 

to u s in the shock-wave case and is equal to -a e in the zero -thickness 
expansion -wave case. 


The following relations result from the equations of motion and 
energy (eqs. (l6) and (17)) for constant potential flow: 


_ „ c a£f , a^f 
^ 1 a-n 3 aq 2 


Sh - ^f + u, £ 

21 Vx2i - 1 ae 


(1 - 1) + 


af 

Hw dq 


iu-j ^ = 0 

1 ^ aqae 


(Bl) 


P-llii C w 


h n ^ /-v 2 

+ (y _ 1)^_ SJ1] 
^Pr ^ 8q^ ^ y^ 2 / 




JL -if if , Sf 1 
21 ■ + Ss j 


h * 

_1 as 


(e - 1) + 


af 


= 0 


(B2) 


Boundary conditions are the same as those given by equations (l8) and the 
coordinates | end q are defined in equations (15); these values are 
valid for both cases under consideration. Note that | is zero on the 
wave in each case. Solutions in the following regions are desired 
(see fig. 9(a)): 

(a) Shock -wave case; region p 


0^|=1 
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(b) Zero -thickness expansion case; region a 


0 < | = 


1 + -Sf s ( 1 + ^) 


Note also that 




and that, in general. 


( h i)„ * Wp 


Direct solutions of equations (Bl) and (B2) appear unlikely. The 
coefficient of the leading partial derivative with respect to | in 
equations (Bl) and (B2) is given by 


(I - 1) + =L 


u 


(e - 1)+ ^ 


and the vanishing of this coefficient introduces singularities into the 
flow field at positions given by 


or 


1 > | £ 1 



0 ^ 


u s t 



u 

— s 


for the flow behind a shock, and 
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or 


0 ^ 


i-Su, 
a c t ~ ^ 


for the case of flow behind a zero-thickness expansion fan. 

Power-series solutions are assumed, as in the expansion-fan case, 
and are given by 

S^U) (B3) 

n=0 


g(i,n) l“gn(n) (A) 

n=0 


Became of the location of the singularities, it is apparent that this 
power series is valid only in the following regions: 

Shock-wave case: 


0^1 < 1 - 


Hi 


Hs 


Expansion-wave cane: 


o £ 6 < l 


For the flow behind the shock, this region is precisely the region (3; 
hence the series solution is sufficient. The region of validity of the 
series solution for the expansion-wave case is only part of the cold- 
gas region a, and it is in this region that the series solution will 
be inadequate . 

Equations (B5) and (b 4) are substituted into equations (Bl) and (B2) 
after expanding and collecting terms, coefficients of like powers of £ 
are equated. The first three equations of each set and the corresponding 
boundary conditions are 
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PiEiVo" + i f o" - | gVo" = o 


^iVi" + (t-| ' ( x - - 1 j^S>* - 0 







- ^-fnf 

2 u* 


l 1 ! 



f 0 (°) = o' 

V(°) = °' 

fo'C*) = 1 


f n (0) = f n '(0) =. f n '(») =0 (e>i) 


filH-i 0 -* 


_ I! vj 

g 0 . (y _ T \ ~1 f H 2 

H Pr V ° 


+ 1®°' " 2 i^ 0 ® 0 ' = 0 


2iHl c w| 


+ 2(r - ^O^l* 


41-^oW- 


1 - t^'k - §!w - ° 






■W 


g + (^" 2 + w) ( -^ 






5 % 

Sl " 5 ^ f2g °’ 


= 0 


(B5a) 

(B5b) 

(B5c) 

(B6a) 

(B6b) 

(BTa) 

(B7b) 

(BTc) 
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s 0 (o) = 

go( w ) = \ 


(B8a) 


gn(0) = 6n(“) = 0 


(a £ l) (B8b) 


Equations (B5a) and (B7a) for Tq and gQ may be integrated numeri- 
cally. Solutions of the remaining equations are 


tqCfl) — ^(^ — * * • “ ® 

= g 2 (n) = . . . = 0 


(B9) 


Thus, both infinite -series representations of the functions reduce to the 
simple relations 


f(S>Tl) = fo(^l) 
g(i,T]) = go(n) 


(bio) 


For constant potential flow, the similar stream and enthalpy functions 
are functions of t) only, a result which could have been obtained directly 
by modifying the similarity transformations (eqs. (13) and (l4)). 

Equations (B5a) and (B7a) and the appropriate boundary conditions may 
be transformed into equations (9) and (10) of reference 9 for both shock 
«nri zero-thickness expansion waves if the special case h^ = h e is consid- 
ered. For this reason, no attempt to evaluate the constant-potential-flow 
solutions is made herein. Singularities apparent in the present solution, 
however, were not evident in reference 9* These singularities represent 
vorticity, generated by the start of motion at t = 0, which diffused 
normal to the wall at an infinite rate and then was washed downstream. By- 
considering both problems as independent of time when treated in a coordi- 
nate system fixed to the wave, reference 9 neglected the start of motion 
and hence, in the shock-tube case, implied a^wave which originated at 
t = propagated through the point x = t_- 0, and continued on. In 
this case singularities were thus not evident. 

Figure (12) shows the location of the singularities. When the flow 
behind the shock wave is considered, only region (3 is of interest; thus. 
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the singularities existing outside this region do not affect the solu- 
tion. Only flow behind the zero thickness expansion wave involves 
singular behavior within the region of interest, and this phenomenon 
should, in principle, be considered in computing shock-tube bo undar y 
layer. The expans ion -wave solution of reference 9 then is act ual ly 
valid only in the region 0 = £ <1. 
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TABLE I.- VELOCITY AMD ENTHALPY FUNCTIONS, n = 0 TEEMS 

[»Pr “ °-T2] 


T 

*0 

f o' 

-p » 
f 0 

so 

_ 11 
So 

0 

0 

0 

1.1284 

0 

0.9575 

.20 

.0213 

.2064 

.9396 

.1775 

.8203 

; . 4 o 

.0802 

.3775 

.7732 

.3290 

.6969 

.60 

.1702 

.5172 

.6284 

.4572 

.5868 

.80 

.2834 

.6301 

.5043- 

.5646 

• 4897 

1.00 

.4208 

.7201 

.3993 

.6539 

.4048 

1.20 

• 5722 

.7910 

.3119 

.7273 

.3314 

i. 4 o 

.7562 

.8459 

.2402 

.7872 

.2687 

l. 6 o 

.9098 

.8880 

.1823 

.8355 

.2157 

1.80 

1.0908 

.9197 

.1364 

.8740 

.1713 

2.00 

I.2772 

.9432 

.1005 

.9045 

.1346 

2.20 

1.4677 

.9604 

.0729 

.9283 

.1047 

2.40 

1.6612 

.9728 

.0521 

.9468 

.0805 

2.60 

1.8567 

.9816 

.0366 

.9609 

.0613 

2.80 

2.0537 

.9878 

.0253 

•9715 

.0461 

3.00 

2.2518 

.9920 

.0172 : 

.9795 

.0343 

3.20 

2.4505 

.9948 

.0115 

.9854 

.0252 

3.40 

2.6497 

.9967 

.0076 

.9897 

.0183 

3.60 

2.8492 

• 9979 

.0049 

.9929 

.0131 

3.80 

3.0489 

.9987 

.0031 

.9951 

.0093 

4.00 

3.2488 

• 9992 

.0020 . 

.9967 

.0065 

4.20 

3.4487 

• 9995 

.0012 

.9978 

.0045 

4.40 

3.6486 

• 9997 

.0007 

.9985 

.0031 

4.60 

3.8486 

.9998 

.0004 

.9990 

.0021 

4.80 

4.0486 

• 9999 . 

.0003 

.9994 

.0014 

5.00 

4.2487 

1.0000 

.0001 

.9996 

.0009 

5.20 

4.4487 

1.0000 

.0001 

.9997 

.0006 

5.40 

4.6488 

1.0000 

.0000 

.9998 

.0004 

5.60 

4.8488 

1.0000 

.0000 

.9999 

.0002 

5.80 

5.0488 

1.0000 

.0000 

.9999 

.0002 

6.00 

5.2489 

1.0000 

.0000 

1.0000 

.0001 

6.20 

5.4489 

1.0000 

.0000 

1.0000 

.0001 

6.40 

5.6490 

1.0000 

.0000 

1.0000 

.0000 

6.60 

5.8490 

1.0000 

.0000 

1.0000 

.0000 

6.80 

6.0491 

1.0000 

.0000 

1.0000 

.0000 

7.00 

6.2491 

1.0000 

.0000 

1.0000 

.0000 




NASA IN 39^3 


39 


TABLE H.- VELOCITy AND ENTHALPY FUNCTIONS, n = 1 TERMS 

[Nj* = 0.72j 7 = 1 . 4 ] 



*1 

f l' 

fl" 

Si 

Si' 

0 

0 

0 

0.7946 

0 

0.1507 

.20 

.0142 

.1336 

.5487 

.0370 

.2084 

AO 

.0505 

.2222 

.3^37 

.0798 

.2125 

.60 

.1006 

.2738 

.1783 

.1198 

.1835 

.80 

.1580 

.2980 

.0502 

.1521 

.1373 

1.00 

.2175 

.2961 

-.0439 

.1743 

.0852 

1.20 

. 275 ^ 

.2805 

-.1080 

.1863 

.0350 

1.40 

.3291 

.2546 

-.1468 

.1888 

-.0086 

l. 6 o 

.3769 

.2231 

-.1652 

.1834 

-.0434 

i. 8 o 

.4182 

.1896 

-.1680 

.1721 

-.0684 

2.00 

.4528 

.1567 

-.1597 

.1567 

-.0843 

2.20 

.4810 

.1262 

-.1443 

.1389 

-.0922 

2.40 

.5035 

.0992 

-.1252 

.1202 

-.0935 

2.60 

.5210 

.0762 

-.1048 

.1018 

-.0899 

2.80 

• 5342 

.0572 

-.0850 

.0845 

-.0830 

3.00 

.5441 

-.0421 

-.0670 

.0688 

-.0740 

3-20 

.5913 

.0303 

-.0514 

.0549 

-.0640 

3.40 

.5564 

.0213 

-.0384 

.0431 

-.0540 

3.60 

.5600 

.0147 

-.0280 

.0333 

-.0445 

3.80 

.5624 

.0100 

-.0200 

.0253 

-.0358 

4.00 

.5841 

.0066 

-.0139 

.0189 

-.0283 

4.20 

.5651 

.0043 

-.0095 

.0139 

-.0219 

4.40 

.5658 

.0027 

-.0063 

.0101 

-.0166 

4.60 

.5663 

.0017 

-.0041 

.0072 

-.0124 

4.80 

.5665 

.0011 

-.0026 

.0050 

-.0091 

5.00 

.5667 

.0006 

-.0016 

.0035 

-.0066 

5.20 

.5668 

.0004 

-.0010 

.0024 

-.0046 

5.40 

.5869 

.0002 

-.0006 

.0016 

-.0032 

5.60 

.5869 

.0001 

-.0003 

.0011 

-.0022 

5.80 

.5869 

.0001 

-.0002 

.0007 

-.0015 

6.00 

.5869 

.0001 

-.0001 

.0004 

-.0010 

6.20 

.5869 

.0000 

-.0001 

.0003 

-.0006 

6.40 

.5869 

.0000 

.0000 

.0002 

-.0004 

6.60 

.5669 

.0000 

.0000 

.0001 

-.0003 

6.80 

.5669 

.0000 

.0000 

.0001 

-.0002 

7.00 

.5670 

.0000 

.0000 

.0000 

-.0001 

7.20 

.5670 

.0000 

.0000 

.0000 

-.0001 

7.40 

.5670 

.0000 

.0000 

.0000 

.0000 

7.60 

.5870 

.0000 

.0000 

.0000 

.0000 

7.80 

.5670 

.0000 

.0000 

.0000 

.0000 

8.00 

.5670 

.0000 

.0000 

.0000 

.0000 




■p- 

o 


oo~j _<| _<i -q ov ov o\ on o\vji vn vji vji vji p- -p- -p" p- -p-vm yi yi yj uj ro ro ro ro ro h h h h h ^ . . . ° 

8'okfe^8§S'Sg8*g > ^fe^ 8§S'S^8g'ogS8§'ofeS8§oS'g 

J 

1 iiiilillllllllll II 1111 

■ , r i ?j i 1 1 i u i i l_i u_i i — j i — it — i u_i l — i u j — i [—j | — • j—j ] — 1 j-j o o o o o o o o o o o o o o o o o o o 

Hi 

ro 

8 *8 8 *8 8 8 8 8 *8 8 8 *8 8 *8 8 8 *8 *8 P p 8888 &S 88 S S888 8 *P 8888 8 

8o88oooooooopHroMvj!oDpo\MO^OMP o\vn o h -p-vo own o o\ own p- -p- 

Hi 

1 1 1 1 1 1 1 1 ■ o 

88888888888888 88PP888 888 

qqqooooooh * h ‘H 1 rv> vh vjh — -3 P vji h -vj vji ro M H O VM VD ON ON -P ON VJI ON H W O O O NO CO 

& 

8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 P P 8 888 888888 88 P 8 8 888888 8 8 

oSoOoBoHHrOrOVHVJl ON VO lO ON O VJI f-* ^ V>4 VO 4=~— 3 -V -F* ON VJI OJhj >3 ON O V>J VJI VJI -VJ IV) 

8* 

1 

1 1 l 1 1 1 1 1 l l J o 

8 8 8 8 8 8 8 8 8 8 8 8 8 P P P 8 8 8888 88888 888 S P g'88P8 fe ! £88 

SoOOOHHMIOyPOvroOUH O -P"~J O H 0-4 O 8—4 ff(5) O H-n P VJI o ro OvVJl VO Ov 

■V 




II 

o 

ro 

v 

n 

H 

. V . 




v_w 

VO 

-p- 

vx 


TABLE III.- VELOCITY AND 




go dT i 


0 




0 


0 



0 


S2 dT l 


.2 

.022 

.4 

.070 

.6 

.149 

.8 

.253 

1.0 

.380 

1.2 

.512 

1.4 

.670 

1.6 

.823 

1.8 

1.000 

2.0 

1.18 

2.2 

1.37 

2.4 

1.54 

2 . 6 ' 

1.74 

2.8 

1.92 

3-0 

2.13 

3.2 

2.32 

3.4 

2.52 

3.6 

2.72 

3-8 

2.92 

4.0 

3.42 

4.2 

3.32 

4.4 

3.52 

4.6 

3.72 

4.8 

3.92 

5.0 

4.12 

5-2 

4.32 

5.4 

4.52 

5.6 

4.72 

5.8 

4.92 

6.0 

5.12 

6.2 

5.32 

6.4 

5.52 

6.6 

5.72 

6.8 

5.92 

7.0 

6.12 


.004 

.015 

.035 

.064 

.095 

.132 

.169 

.208 

.244 

• 275 
.304 
.330 
.552 

• 572 
.587 
.598 
.408 
.416 
.422 

.426 

.429 

.431 

.432 

.434 

.435 

.436 

.436 

.437 

• 437 

• 437 
.437 
.437 
.437 
.437 
.437 


.000 

,003 

.011 

.025 

.041 

.057 

.071 

.082 

.089 

.091 

.088 

.084 

.075 

.065 

.054 

.042 

.031 

.022 

.012 

.004 

-.002 

-.008 

-.013 

-.016 

-.019 

-.021 

-.023 

-.024 

-.024 

-.025 

-.025 

-.026 

-.026 

-.026 

-.026 
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(a) Wave diagram. 



Figure 2.- Shock- tube wave diagram and boundary layer for strong expan- 
sion fan. 
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nth-order velocity function, f n ‘ ( -17 ) 


(a) Velocity functions. 
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nth-order enthalpy f u n cti on, g n ( 77) 


(b) Enthalpy functions. 

Figure 3 .- Velocity and enthalpy function profiles of the nth order, n = 0 , 

1 , and 2. 
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Figure 4 


Enthalpy 


difference 


ratio 



(b) Enthalpy difference profiles. 

- Velocity and enthalpy difference profiles for various values of 
the conical parameter | . 
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MCA TN 39^5 



(a) Skin friction. 

Figure 5.- Skin-friction and he at -transfer functions near the leading 

edge of the expansion fan. 
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Figure 6.- Variation of compressible similarity parameter with incompres 

sible parameter." 
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MCA TN 39*4-3 



Figure 7 «- Concluded. 


're e - stream stagnation enthalpy difference, 






Conical coordinate, £ 


(a) Skin friction. 

Figure ID.- Skin-friction and heat-transfer functions in the expansion fan compared with those 
for a zero-thickness expansion wave. The asterisks represent the appropriate positions of 









Heat-tr a n sfer coefficient, q Skin-friction coefficient, 
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(a) Skin friction. 



Conical coordinate at expansion-fan trailing edge , £f e 


(b) Heat transfer. 

Figure 11.- Skin-friction and heat-transfer coefficients for expansion 
fan and for zero-thickness expansion wave. 









